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Abstract 

Performance  characteristics  of  micro,  single-chamber,  intermediate  temperature  solid  oxide  fuel  cell  (IT-SOFC)  has  been  successfully 
evaluated  by  computational  simulation  based  on  macro  modelling.  The  system  has  composite  electrodes  that  are  alternately  placed  on  the 
ceria-based  electrolyte.  The  polarization  curve  is  obtained  through  the  series  of  computational  processes  for  structural  parameters  such  as 
distance  between  the  anode  and  the  cathode  and  the  thickness  of  the  electrolyte.  Also,  the  effect  of  the  inflow  direction  of  fuel-air  mixture  on 
the  performance  of  the  fuel  cell  is  quantitatively  investigated. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  solid  oxide  fuel  cell  (SOFC)  is  highly  efficient  energy- 
conversion  system  that  converts  chemical  energy  to  electrical 
energy.  It  can  utilize  hydrocarbons  such  as  CH4,  C2H6  and 
C3H8  and  CO  besides  H2  as  a  fuel  source.  No  additional 
fuel-reformer  system  is  required  because  of  the  high  oper¬ 
ation  temperature.  The  SOFC  can  be  applied  to  versatile 
power  generation  systems,  such  as  stationary  power  supplies, 
auxiliary  power  sources  for  automobiles  and  military  appli¬ 
cations.  Research  has  been  mainly  focused  on  the  fabrication 
of  various  forms  of  SOFC,  e.g.,  tubular  and  planar  types,  for 
commercial  applications.  Experimental  attempts  have  also 
been  made,  [1,2]  to  develop  develop  micro  SOFCs  (jjlSOFCs) 
as  power  sources  of  portable  equipment  and  MEMS  appli¬ 
cation  [1,2].  To  achieve  stability  and  economical  efficiency, 
intermediate  temperature  SOFCs  (IT-SOFCs)  with  ceria- 
based  electrolyte  have  been  proposed  [3-6] .  Materials  such 
as  stainless  steel  can  be  employed  as  a  bipolar  plate  in 
IT-SOFC  systems  and,  consequently,  low  processing  cost 
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can  be  attained.  Gas  sealing  has  been  a  vital  issue  and  the  gas¬ 
sealing  amorphous  materials  must  be  selected  on  the  basis  of 
the  degree  of  thermal  expansion  matching  between  compo¬ 
nents  and  mechanical  strength.  In  this  respect,  it  is  interesting 
to  note  that  gas- sealing  problems  can  be  largely  mitigated  by 
using  a  single-chamber  SOFC  system  since  no  separation 
between  fuel  and  air  is  required  [7-12]. 

In  a  SOFC  system,  computer  simulation  based  on  the¬ 
oretical  modelling  is  known  to  be  a  very  efficient  method 
to  predict  materials  properties,  reaction  kinetics  and  unit 
cell  performance.  Theoretical  modelling  of  SOFC  can  be 
divided  into  three  parts.  First,  basic  properties  such  as  the 
ionic  conductivity  of  electrolyte  materials  can  be  predicted 
by  molecular  dynamics  or  ab  initio  approaches  [13-17].  Sec¬ 
ond,  the  reaction  kinetics  at  the  three-phase  boundary  (TPB) 
at  the  electrode  can  be  simulated  by  micro  modelling  [  1 8-2 1  ] . 
Third,  a  detailed  understanding  of  the  fluid  dynamics,  mass 
transport  and  heat  conduction/convection  in  the  SOFC  system 
can  be  obtained  through  macro  modelling.  The  overpoten¬ 
tials  of  individual  SOFC  cells  can  be  successfully  predicted 
through  a  series  of  macro  modelling  procedures  [22-25]. 

In  this  study,  fluid  motion,  diffusion  of  the  fuel-air  mix¬ 
ture  and  oxygen  ion  conduction  is  considered  in  the  macro 
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Nomenclature 

Ci 

concentration  of  species  i  (kgm-3) 

Di 

diffusivity  of  species  i  (m2  s-1) 

F 

Faraday  constant  (C  mol-1) 

channel 

concentration  of  species  i  at  the  gas  channel 
(mol  m-3) 

h 

current  density  (Am-2) 

k 

exchange  current  density  (Am-2) 

mTPB 

concentration  of  species  i  at  three  phase  bound¬ 
ary  (molm-3) 

Mi 

molar  weight  of  species  i  (kg  mol-1) 

N0.i 

inward  mass  flux  of  species  i  (kgm-2  s-1) 

P 

pressure  in  channel  (Pa) 

P2 

pressure  in  electrodes  (Pa) 

Pn-S 

Pressure  obtained  from  calculating  the 
Navier-Stokes  equation  (Pa) 

R 

gas  constant  (J  mol-1  K-1) 

Rt 

volumetric  reaction  rate  of  species  i 
(kgm-3  s-1) 

T 

temperature  (K) 

u 

velocity  vector  (ms-1) 

U ,  V,  w 

x-,  y-,  z-components  of  velocity  in  channel 
(ms-1) 

U2,V2,W2  x-,y-,  z-components  of  velocity  in  electrodes 
(ms-1) 

v; 

stoichiometric  coefficient  of  species  i 

v 

voltage  (V) 

Kact 

activation  loss  (V) 

I^con,a 

concentration  loss  in  anode  (V) 

^con,c 

concentration  loss  in  cathode  (V) 

Greek  letters 

Q'a 

charge  transfer  coefficient  at  anode 

OL  c 

charge  transfer  coefficient  at  cathode 

0 

dynamic  viscosity  (Pa  s) 

K 

permeability  (m2) 

P 

density  (kg  m-3) 

a 

conductivity  (S  m-1) 

modelling  scheme  by  using  the  three-dimensional  finite 
element  method  (FEM).  Next,  ohmic  loss,  activation  loss 
and  concentration  loss  are  obtained  through  a  series  of  macro 
modelling  procedures  and,  consequently,  the  performance 
characteristics  of  a  micro  single-chamber  IT-SOFC  are 
quantitatively  investigated  for  various  structural  parameters. 


2.  Computational  process 

For  the  simulation,  a  Hibino  B-type  single-chamber  SOFC 
[9,12]  was  employed  (Fig.  1).  The  dimensions  of  the  model 
for  computer  simulation  are  listed  in  Table  1 . 
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Fig.  1 .  Schematic  diagram  of  calculation  model  based  on  Hibino  B-type  fuel 
cell  [9,12]:  (a)  side  view  and  (b)  top  view. 

A  fully  dense  and  non-porous  electrolyte  was  assumed 
and  no  penetration  of  gas  species  into  the  electrolyte  was 
allowed.  A  NiO-Ceo.9Gdo.10 1.95  (Ni-GDC)  anode  and  a 
Smo.5Sro.5Co03-10%Sm2Ceo.801.9  (SSC-10SDC)  cathode 
with  uniform  pores  of  100  nm  diameter  were  designed  to 
exhibit  electronic  and  ionic  conducting  properties.  It  was 
also  assumed  that  the  oxidation  of  fuel  in  the  gas  mixture 
and  the  reduction  of  oxygen  occurred  only  at  the  anode  and 
at  the  cathode,  respectively.  The  TPB,  where  the  electrode 
reaction  occurs,  was  assumed  to  be  the  interface  between  the 


Table  1 


Physical  dimensions  of  micro  single-chamber  SOFC 


Symbol 

SOFC  component 

Dimension  (pan) 

E_L 

Electrolyte  length 

1000 

E_W 

Electrolyte  width 

160 

E_T 

Electrolyte  thickness 

20 

A_L 

Anode  length 

900 

A_W 

Anode  width 

20 

A_T 

Anode  thickness 

20 

C_L 

Cathode  length 

900 

C_W 

Cathode  width 

20 

C_T 

Cathode  thickness 

20 

G_H 

Height  of  hypothetical  gas  channel 

200 
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electrolyte  and  electrode,  because  most  electrode  reactions 
have  been  found  [24]  to  take  place  within  a  50  p,m  range  in 
the  interface  region  of  the  electrolyte  and  the  electrode.  In 
the  present  modelling,  however,  the  TPB  was  assumed  to 
be  the  whole  electrode  volume  because  the  thickness  of  the 
electrodes  was  less  than  50  pan. 

For  the  simplicity,  only  hydrogen  fuel  was  considered.  The 
respective  electrode  chemical  reactions  can  be  represented 
by: 

Anode  :  2H2  +  02  =  2H20  +  4e"  (1) 

Cathode  :  02  +  4e~  =  202"  (2) 

The  injected  fuel-air  mixture  was  composed  of  hydrogen, 
nitrogen,  oxygen  and  a  small  amount  of  water  vapour  and 
the  temperature  of  the  unit  cell  system  was  isothermally 
fixed  at  773  K.  A  series  of  finite  element  calculations  to 
find  solutions  for  ionic  conduction,  fluid  motion,  diffusion 
and  convection  partial  differential  equations  were  carried 
by  FEMLAB  v3.1.  The  ohmic,  activation  and  concentration 
losses  were  obtained  from  calculation  of  the  voltage  drop  due 
to  ionic  conduction,  finding  the  roots  of  the  Butler- Volmer 
equation  and  determining  the  concentration  distribution  of 
each  species  resulting  from  the  diffusion  and  convection 
in  the  anode,  cathode  and  gas  channel.  The  detailed  pro¬ 
cedure  for  modelling  performance  characteristics  from  loss 
calculations  is  depicted  in  Fig.  2.  Finally,  the  theoretical  open- 
circuit  voltage  (OCV)  was  obtained  from  the  Nernst  equation 
[26]. 


2.7.  Ohmic  loss 

The  ohmic  loss  is  the  voltage  drop  caused  by  the  resistance 
of  migrating  oxygen  ions  and  electrons  in  the  electrodes.  In 
this  work,  only  ohmic  loss  from  ionic  conduction  has  been 
considered  since  that  originating  from  electronic  conduction 
is  negligible.  Using  reported  values  of  the  conductivity  (a)  of 
oxygen  ions  in  the  electrolyte  [27]  and  the  electrodes  [27,28] 
at  773  K,  Eq.  (3)  was  solved,  numerically  and  the  voltage 
obtained. 

-V  •  (crVV)  =  0  (3) 

The  boundary  conditions  applied  to  the  calculation  are  given 
by  Eqs.  (4)— (6)  for  inward  current  density,  outward  current 
density  and  electrical  insulation: 


-n  •  (-crVV)  =  /e 

(4) 

— n  •  (—crVV)  =  -ie 

(5) 

— n  •  (—crVV)  =  0 

(6) 

2.2.  Activation  loss 

The  activation  loss  is  due  to  overcoming  the  energy  barrier 
when  electrons  released  during  the  electrochemical  reac¬ 
tion.  Therefore,  the  activation  loss  depends  on  temperature, 
fuel  species,  electrode  materials,  oxygen  partial  pressure, 
etc.  There  are  many  different  opinions  of  the  nature  oxida¬ 
tion  and  reduction  mechanisms  in  SOFC  electrodes,  hence 
precise  elucidation  of  multi-step  electrochemical  reactions 
are  still  remained  unsolved  [22].  The  activation  loss  of 
IT-SOFC,  however,  becomes  more  critical  compared  with 


Fig.  2.  Flow  chart  of  performance  characteristics  analysis  based  on  macro  modelling. 
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high  temperature-SOFCs.  The  Butler- Volmer  Eqs.  (7a)  and 
(7b)  represents  the  relation  between  activation  loss  and  cur¬ 
rent  density: 


4 


exp 


—  exp 


(7a) 


+  0ic  —  1 


(7b) 


Exchange  current  density  (z'o),  charge  transfer  coefficient  (aa 
and  ac)  were  referred  to  results  of  Zhu  and  Kee  [22].  The 
numerical  solution  for  the  activation  loss  ( Vact)  was  obtained 
by  the  Newton-Raphson  method  [29]. 


2.3.  Concentration  loss 


applied  to  other  boundaries  of  electrode.  The  resulting  pres¬ 
sure  values  from  fluid  calculations  at  the  gas  channel  were 
stored  and  used  as  a  boundary  value  for  the  porous  electrodes 
Eq.  (14). 


u  (  KVP  )  -0 

(13) 

V  v  ) 

p  =  Pn-S 

(14) 

Finally,  from  the  velocity  distribution,  the  steady-state  diffu¬ 
sion  and  convection  equation,  Eq.  (15)  was  solved  to  obtain 
the  concentration  distribution  of  reactants  and  products  in  the 
gas  channel  and  the  electrodes. 

V  •  (- Di Vet  +  Cjii)  =  Ri  (15) 


Differences  in  the  concentration  of  each  species  between 
the  gas  channel  and  the  TPB  lead  to  a  concentration  loss.  In 
previous  studies  [24],  the  gas  velocity  in  the  porous  electrodes 
was  set  to  zero.  In  this  work,  however,  convective  diffu¬ 
sion  was  determined  by  calculating  the  velocity  of  gas  fluid 
inside  both  electrodes.  The  inlet  gas  mixture  was  composed 
of  H2:N2:02:H20  =  0.5:0.38:0. 105:0.0105  in  mole  fractions 
and  assumed  to  obey  the  ideal  gas  law  and  amount  of 
1 80  seem  (H2 : 90  seem,  air:  90  seem)  was  injected  into  the  gas 
channel.  To  calculate  the  velocity  distribution  in  the  gas  chan¬ 
nel,  the  steady-state  incompressible  Navier-Stokes  equation, 
Eq.  (8a)  with  Eq.  (8b)  was  solved  numerically: 

p(u  •  Vu)  —  V  •  77(Vu  +  (Vu)T)  +  V p  =  0  (8a) 

V  •  u  =  0  (8b) 

The  density  (p)  was  calculated  from  the  ideal  gas  law  and  the 
dynamic  viscosity  ( r) )  of  the  fuel-air  mixture  was  the  molar 
mean  value  obtained  from  the  product  of  individual  viscosity 
and  mole  fraction.  For  setting-up  boundary  conditions,  an 
initial  gas  velocity  (uq)  of  42.47  m  s-1  and  a  pressure  (Pq)  of 
101325  Pa  were  applied  to  the  gas  inlet  and  outlet  boundaries 
Eqs.  (9)  and  (10).  At  the  surface  of  electrolyte  and  electrode, 
no  slip  conditions  were  employed  Eq.  (11). 

u  n  =  uo  (9) 

n  •  (—V  •  ?7(Vu  +  (Vu)T)  +  VP)  =  — nPo  (10) 

u  =  0  (11) 

To  obtain  the  velocity  distribution  in  the  porous  elec¬ 

trodes,  Darcy’s  law  equation  (Eq.  (12))  was  solved  with 
the  calculated  pressure  values  that  resulted  from  solving  the 
Navier-Stokes  equation. 

V-(  —  VPj=0  (12) 

The  permeability  (/c)  of  the  anode  and  the  cathode  was 
assumed  to  be  1.0  x  10“ 12  m2.  The  normal  flux  of  the  mix¬ 
ture  gas  was  set  to  zero  Eq.  (13)  at  the  interface  between  the 
electrolyte  and  the  electrode  and  the  pressure  condition  was 


The  value  of  the  diffusivity  (Di)  of  each  species  in  the  chan¬ 
nel  and  porous  media  is  given  in  Table  2.  The  diffusivity 
was  calculated  from  the  Hirschfelder  model  based  on  the 
Lennard-Jones  potential  [30].  In  particular,  when  the  diffu¬ 
sivity  in  the  porous  electrode  was  calculated,  Knusen  dif¬ 
fusion  as  well  as  the  Bruggman  model  was  considered  to 
represent  the  actual  porous  media,  which  consisted  of  ran¬ 
domly  distributed  tortuous  pores.  The  volumetric  reaction 
rate  (Ri)  was  defined  as  a  mass  quantity  of  reactants  and 
products  per  unit  volume  per  second,  i.e., 

Ri  =  zb -  x  area  of  reaction  site(ARS) 

4  F 

1 

x -  x  Mi  (16) 

electrode  volume 

where  Ri  is  positive  if  species  i  is  the  product  (H2O)  and  is 
negative  if  species  i  is  the  reactant  (H2,  O2);  v*  is  the  stoichio¬ 
metric  coefficient  of  reactions,  i.e.,  vh2  =  vh2o  =  2,  vo2  = 
1 .  The  area  of  the  reaction  site  (ARS)  can  be  calculated  from 
the  permeability. 

The  mass  flux  obtained  from  the  density  and  the  influx 
volume  of  the  gas  mixture  was  applied  to  the  inlet  boundary 
as  a  boundary  condition  Eq.  (17).  At  the  interface  between 
the  electrolyte  and  the  electrodes,  migration  of  each  molecule 
across  the  interface  should  not  be  allowed  Eq.  (18).  A  con¬ 
vective  flow  condition  was  applied  to  the  outlet  boundary  Eq. 

(19). 


n  •  (-D,Vc,  +  c,-u)  =  N0J 

(17) 

n  •  (— A'Vc,-  +  c,-u)  =  0 

(18) 

n  •  {-DiVa)  =  C[U  ■  n 

(19) 

Table  2 


Diffusivity  of  each  species  at  773  K,  1  atm  (cm2s  !) 


Channel 

Porous  electrode 

H2 

2.813 

0.252 

n2 

2.366 

0.081 

o2 

0.559 

0.059 

h2o 

0.757 

0.079 
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From  the  calculated  concentration  distribution,  the  con¬ 
centration  loss  was  calculated  through  Eqs.  (20)  and  (21) 
[22]. 


RT  (  [H2] 

W  =  —  I  2  In 


channel 


4  F 


[H2] 


TPB 


—  2  In 


[H20] 


channel 


[H20] 


TPB 


(20) 


V  _  RT  I ,  [02] 

Fcon,c  —  '  In 


channel 


4  F 


[02] 


TPB 


(21) 


3.  Results  and  discussion 

The  calculated  loss  corresponding  to  variation  in  the  cur¬ 
rent  density  is  shown  in  Fig.  3a.  As  with  other  types  of  SOFC, 
the  major  loss  originates  from  the  ohmic  loss  caused  by  the 
diffusion  of  oxygen  ions.  On  the  other  hand,  the  concentration 
loss  is  minimal  due  to  the  reduced  concentration  difference 
between  the  inside  of  the  electrodes  and  the  gas  channel  since 
whole  volume  of  the  electrodes  served  as  reaction  sites.  The 
corresponding  polarization  curve,  which  is  directly  related 
to  the  performance  of  the  SOFC,  was  obtained  by  subtract¬ 
ing  all  the  calculated  losses  from  the  theoretical  open-circuit 
voltage  (1.133  V)  (see  Fig.  3b). 
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Fig.  3.  (a)  Predicted  losses  as  current  density  variation,  (b)  Predicted  polar¬ 
ization  curve  the  distance  between  electrodes  and  the  thickness  of  electrolyte 
is  20  and  20  pan,  respectively. 
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Fig.  4.  Polarization  curve  corresponding  to  variation  in  the  distance  between 
two  electrodes. 


Next,  the  effect  of  the  distance  between  the  anode  and  the 
cathode  on  the  polarization  curve  was  quantitatively  investi¬ 
gated  (Fig.  4).  On  decreasing  this  distance,  the  power  density 
of  cell  increases  due  to  a  reduction  in  ohmic  loss.  This  agrees 
well  with  empirical  expectation  and  is  consistent  with  exper¬ 
imental  reports  [9,12]. 

The  polarization  curve  according  to  variation  in  the  thick¬ 
ness  of  the  electrolyte  is  presented  in  Fig.  5.  In  contrast  to 
a  Hibino  A-type  single-chamber  SOFC  [9,12],  the  single¬ 
chamber  SOFC  adapted  in  this  study  (B-type)  gives  a  higher 
power  density  with  a  thicker  electrolyte.  This  is  due  to  the 
ohmic  loss  behaviour  according  to  variation  in  the  thickness 
of  the  electrolyte  (see  Fig.  6).  The  ohmic  loss  decreases  with 
increase  in  the  thickness  and  converges  at  thicknesses  above 
80  pim.  It  can,  therefore,  be  reasonably  concluded  that  in 
the  manufacturing  step  of  a  B-type  single-chamber  SOFC, 
a  thicker  electrolyte  can  be  employed  without  deteriorating 
the  performance  characteristics,  and  no  additional  support 
materials  for  the  required  mechanical  strength  are  required. 

The  effect  of  variation  in  the  direction  of  the  inflow  gas 
mixture  on  the  polarization  curve  is  demonstrated  in  Fig.  7. 
The  highest  power  density,  is  obtained  when  the  gas  mix¬ 
ture  is  injected  into  the  vicinity  of  the  cathode  in  a  direction 


Fig.  5.  Polarization  curve  for  variation  of  the  electrolyte  thickness. 
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Fig.  6.  Ohmic  loss  variation  with  electrolyte  thickness. 


Fig.  7.  Polarization  curve  for  various  inflow  directions  of  fuel-air  gas  mix¬ 
ture. 

perpendicular  to  the  electrode  axis.  This  observation  can  be 
understood  by  considering  the  diffusivity.  The  diffusivity  of 
O2  in  the  cathode  is  much  lower  than  that  of  H2  and  H2O 
in  the  anode,  and,  therefore,  the  migration  of  oxygen  in  the 
cathode  is  significantly  retarded.  By  contrast,  when  inflow 
gas  is  injected  into  the  cathode  side,  the  velocity  of  gas  in 
the  cathode  is  higher  than  that  in  the  anode.  This  increment 
of  fluid  velocity  enhances  the  convective  diffusion  of  oxy¬ 
gen  in  the  cathode  and,  eventually,  lowers  the  concentration 
loss.  On  the  other  hand,  when  the  inflow  gas  direction  is  par¬ 
allel  to  the  electrode  axis,  the  power  density  is  predicted  to 
be  significantly  lowered  because  the  concentration  difference 
of  the  gas  species  between  the  inlet  and  outlet  regions  of  the 
SOFC  system  is  extensive.  This  originates  from  the  fact  that 
the  length  of  both  electrodes  is  45  times  longer  than  the  width 
and  height. 

4.  Summary  and  conclusions 

Computer  simulation  based  on  the  macro  modelling 
of  multi-physical  phenomena,  e.g.,  ionic  conduction,  fluid 
dynamics  and  gas  diffusion,  has  been  successfully  performed. 


Consequently,  ohmic  loss,  polarization  loss  and  concen¬ 
tration  loss  are  quantitatively  derived  from  the  simulation 
results.  The  performance  of  a  micro  scale  single-chamber 
SOFC  system  is  found  to  be  largely  dependent  on  the  struc¬ 
tural  parameters,  the  distance  between  the  electrodes,  the 
thickness  of  the  electrolyte  and  the  inflow  direction  of  the 
gas  mixture.  In  general,  an  improvement  in  performance  can 
be  obtained  with  a  narrow  interdistance  between  the  elec¬ 
trodes  and  with  a  thicker  electrolyte.  When  the  inflow  direc¬ 
tion  is  perpendicular  to  the  electrode  axis,  the  specific  micro 
single-chamber  SOFC  cell  delivers  the  best  performance.  It 
is  expected  that  the  computational  scheme  developed  in  this 
study,  can  serve  as  an  important  guideline  for  manufacturing 
the  micro  single-chamber  IT-SOFC  system  and  for  optimiz¬ 
ing  the  performance  characteristics. 
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